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- - - We present a detailed numerical study of the elementary condensation events (avalanches) asso- 

ciated to the adsorption of ''He in silica aerogels. We use a coarse-grained lattice-gas description 
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and determine the nonequilibrium behavior of the adsorbed gas within a local mean-field analysis, 
neglecting thermal fluctuations and activated processes. We investigate the statistical properties of 
the avalanches, such as their number, size and shape along the adsorption isotherms as a function 
of gel porosity, temperature, and chemical potential. Our calculations predict the existence of a line 

Oof critical points in the temperature-porosity diagram where the avalanche size distribution displays 
I a power-law behavior and the adsorption isotherms have a universal scaling form. The estimated 

, critical exponents seem compatible with those of the field-driven Random Field Ising Model at zero 

temperature. 

Ch ' PACS numbers: 64.70.Fx, 64.60.Ht, 75.60.Ej, 67.70. +n 
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I. INTRODUCTION 



■ This work is part of a series of papers devoted to the study of gas adsorption and capillary condensation in silica 
aerogels^, Q- Our goal is to provide a theoretical interpretation of ''He adsorption experiments that have been 
, performed over the last 15 years in several low-temperature physics groups 0, 0,1^ 0]. More generally, we wish to 
^ ' propose a theoretical framework for describing the nonequilibrium (hysteretic) behavior of fluids in disordered porous 
solids0- It is worth noting that the experimental studies of helium in aerogel were originally aimed at understanding 
^ ■ the effect of quenched random impurities on equilibrium critical phenomena and phase separation. In particular, the 
[ study of the liquid- vapor transition of ^Hc followed from the suggestion that the critical behavior of mixtures or pure 
fluids confined in porous media or in contact with a gel network could be interpreted in terms of the random field 
Ising model (RFIM) ,8] . In this perspective, the use of aerogel was motivated by its very tenuous structure (which 
J — minimizes the effects of confinement) and by the fact that the porosity (and thus the amount of disorder) can be 
\^ varied in a controlled way. However, in spite of the early claim that mesoscopic equilibrium phase separation of "^He 
fS| could be observed in a 95% porosity aerogelQ, leading to an estimation of the order parameter critical exponent /3, 
OO it now seems more and more questionable that true equilibrium behavior is observed in experiments. In particular, 
recent measurements of sorption isotherms using a high-precision capacitive technique^ show that hysteresis between 
filling and emptying persists to temperatures very close to the critical temperature of pure helium (Tc = 5.195/^), 
and gradual adsorption is always observed above T = A.SSK, instead of the sharp vertical step expected for two-phase 
"j^ ■ coexistence. Similar observations have been made by another group 6], using a low- frequency mechanical pendulum 
technique and, more recently, coupling pressure and optical measurements. These results are reminiscent of the 
, behavior of fluids confined in dense porous glasses which is usually interpreted in terms of capillary condensation. 
'"^ ' Although this moves away the exciting perspective of studying the putative equilibrium vapor-liquid critical point 
^ of "'He inside aerogel, the remarkable changes in the hysteretic behavior of the adsorbed fluid with temperature and 
O porosity remain to be explained. In particular, the hysteresis loops in very light aerogels have a rectangular (instead 
of triangular) shape at very low temperature (for instance, at T = 2.3'iK in 98% porosity aerogel|^) with a very steep 
adsorption branch that suggests the presence of a genuine, but nonequilibrium, first-order phase transition. This is 
in agreement with our theoretical description'!', which predicts the change from a continuous to a discontinuous 
adsorption isotherm as one decreases the temperature (at constant porosity) or increases the porosity (at constant 
temperature), the jump in the fluid density corresponding to a macroscopic avalanche in the system. (Note that 
we only focus here on adsorption. The behavior on desorption is also interesting but is related to other physical 
processes as discussed in Refs.P, 0].) This scenario, which resembles the one predicted for the field-driven RFIM 
at zero temperature^, and observed in some disordered magnetic systems^, 0|, suggests the existence of a 
line of disorder-induced critical points in the temperature-porosity diagram with an associated scaling behavior. In 
the present work we further elaborate on this issue by presenting new results concerning the avalanche properties 
as a function of temperature and porosity and by studying the scaling behavior of the adsorption isotherms in the 
vicinity of the nonequilibrium critical points. We are also interested in characterizing the avalanches geometrically 
and studying their possible relation with the aerogel structure. 

The paper is arranged as follows. In section 2 we review the model and the theory. In section 3 we present the 
results for the statistical properties of avalanches. In section 4 we discuss the scaling behavior of the adsorption 
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isotherms and present the phase diagram of the model in the temperature-porosity plane. We conclude, in section 5, 
with a discussion of the relevance of our results to experiments. 



II. MODEL AND THEORY 

Since both the model and the theory were described in detail in ealier work (see in particular Ref.'l'l), we will 
restrict the presentation to some key aspects that are relevant to our present goal. 

The gel-fluid system is modelled by a coarse-grained lattice-gas where each cell of the size of a silica particle 
(typically a few nanometersfH) is occupied either by the solid or the fluid. The distribution of silica particles on 
the lattice is obtained by a diffusion-limited cluster-cluster aggregation (DLCA) algorithm 13J, so as to model the 
intricate microstructure of actual base-catalized silica aerogels. On the other hand, fluid particles can equilibrate with 
a reservoir that fixes their chemical potential /i and the temperature T . Multiple occupancy of a site is forbidden and 
only nearest-neighbors (nn) attractive interactions are considered. For reasons that are explained in Ref.0, we use a 
body-centered cubic (bcc) lattice and periodic boundary conditions (TV = 2L^ where L is the linear size). 

For a given configuration of the solid matrix, fluid configurations are obtained by minimizing a grand-potential 
functional in the local mean field approximation 0, 

^[{Pi}] = kBT^[p^ Inp, + {^, - p,) \n{^, ~ p,)] - wjf ^ piPj - w^f ^ [pt{l ~ r]j) + Pj{l - 77^] - p.'^Pi (f) 

i <ij> <ij> i 

where pi {i = f ...A^) is the thermally averaged fluid density at site i and rji = 0, 1 is the quenched occupancy variable 
for the gel particles {rji = if the site is occupied by the gel and f otherwise), w/f and Wsf denote the fluid-fluid 
and solid-fluid interactions, respectively, and the double summations run over all distinct pairs of nn sites. The gel 
porosity is given by = (l/iV) rji and the ratio y — Wsj/wff controls the wettability of the solid surface. 

f2[{pi}] thus deflnes a grand-potential landscape whose extrema are obtained by solving the corresponding Euler- 
Lagrange equations 6fl/Spi = 0. At flxed T and p, and for a given realization of the solid, this yields a set of N 
coupled nonlinear equations 



where (3 — l/{kBT) and the sum runs over the nearest neighbors of site i. A crucial feature of these equations is 
that they may have a large number of solutions at low temperature (and in a certain range of p) as a consequence 
of the randomness of the solid matrix Ti]. We believe that this feature is the key for understanding the physics of 
adsorption in disordered porous solids. Especially important are the local minima of r2[{/9i}] which are obtained by 
solving the set of Eqs. (2) iterativelyfl^ and which correspond to metastable states (to be distinguished from the 
global minimum of ri[{/3i}] that deflnes the true equilibrium state). Note that despite its mean- field character (by 
which thermal fluctuations are only crudely described), the present approach fully accounts for the disorder- induced 
fluctuations. In contrast with the standard global mean-field approximation, our treatment still allows to describe 
strongly inhomogenous configurations and to predict non-classical critical behavior. 

Computing the metastable states is however not sufficient for describing the nonequilibrium behavior of the system, 
and one also needs to specify a dynamics. Our theoretical descriptionjJI^, Q is based on the use of a non-ergodic, 
deterministic, zero-temperature-\ike dynamics which neglects all thermally activated processes that would allow the 
system to cross the (grand-potential) barriers: a change in the fluid configurations is thus only due to a change in 
the landscape that comes itself from a change in the temperature or from a variation of the external field (here the 
chemical potential p). As p varies (for instance, when the pressure in the reservoir is slowly increased at constant 
T), the system either follows the minimum in which it was trapped as this minimum deforms gradually, or it falls 
instantaneoulsy into another minimum when the former reaches its stability limit. This latter move is a discontinuous 
and irreversible process, an avalanche, which is at the origin of the history-dependent behavior of the system. The 
avalanche corresponds to some collective condensation event inside the gel network which manifests itself by a jump 
in the adsorption isotherm (whereas the gradual deformation of a minimum only corresponds to a slight swelling of 
the liquid domains). The size of the avalanche is thus given by the size of the jump, i.e. the discontinuous change in 
the average fluid density Pf — i^/N)^^ Pi- As was shown qualitatively in Ref.Q, these condensation events occur 
at different length scales inside the gel, depending on chemical potential, temperature, and porosity. In particular, 
when (f) is large enough, for instance = 95%, there exists a critical point separating a high-temperature regime 
where all condensation events are of microscopic or mesoscopic size (and the adsorption isotherm is smooth in the 
thermodynamic limit), and a low-temperature regime where a macroscopic fraction of the gas condenses, which yields 
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a finite jump in the isotherm. It is one goal of the present work to put these observations on a more quantitative 
basis. 

We shall discuss below in more detail the relevance of this theoretical description to the actual experimental 
situation. Let us just recall that this picture is expected to be valid when there are widely separated time scales 
in the problem. Specifically, the observation time must be much larger than the local equilibration time (i.e. the 
time needed to relax from a marginally stable state to a nearby metastable state) but much smaller than the time 
associated to thermal activation (i.e. the time that it takes to escape from a local minimum via a nearby saddle 
point) Jj^]. This requires that energy barriers are much larger than temperature, which is indeed the case in many 
experimental situations dealing with random systems. It may also be noted that our definition of avalanches generalizes 
the one used at T = (see e.g. Ref.[l7|l. We indeed introduce an explicit temperature-dependence in the problem 
by allowing the (mean-field) free-energy landscape itself to change with T (for a similar approach in the context of 
random magnets, see e.g. Ref. |l8j). As shown in our previous works, this dependence explains quite convincingly 
the variations observed in the hysteretic behavior of gases adsorbed in disordered solids at low enough temperature. 

III. AVALANCHE PROPERTIES 

As discussed in the preceding papers of this series[ll0|, the number, size and shape of the condensation events 
along an adsorption isotherm depend on porosity and temperature. They also depend on the value of the chemical 
potential at which they occur. The random character of these events (each finite-size sample has a different isotherm) 
requires a statistical study of their properties, an issue that has been investigated in great detail in the context of the 
field-driven RFIM[l3,|li. 

A. Numerical procedure 

The first problem is to identify each avalanche properly. At T = 0, when the microscopic variables can only take 
discrete values (such as the Ising spins in the RFIM) , this is an easy task because the energy landscape only changes 
through the disappearance of minima, which induces the avalanches. One can then vary the external field until a 
spin becomes unstable (this spin is unique when the random field distribution is continuous) and then count the 
number of spins that flip during the avalanche to determine its size. At finite temperature, on the other hand, one 
has to deal with continuous variables (the local fluid densities pi in the present case), and since Sfi, the increment 
in the chemical potential along the adsorption isotherm, is always finite, one has to distinguish between the gradual 
deformation of a minimum and the jump associated to the passage to another local minimum, both leading to a 
change in the total density. Moreover, if 6fj, is not small enough, several avalanches may overlap, which biases the 
statistics of their number or size distribution. Although, in principle, one could precisely locate the beginning of each 
avalanche by following the evolution of the eigenvalues of the Hessian matrix Hij = d^ft/dpidpj with p, or T (an 
avalanche starts when a local spinodal corresponding to the vanishing of a certain eigenvalue is reached[23), this is a 
daunting numerical task for a large system (with N ~ 10^) and it is more advisable to use another procedure. 

Let us first note that the determination of the location and size of the avalanches depends on the precision in the 
computation of the metastable states. Since these latter are obtained by solving the set of Eqs. (2) iteratively, it must 
be checked that the results are not significantly modified when the convergence criterium of the iterative procedure 
is changed. All the results presented here were computed with a precision e = 10~^ (i.e., the calculation was stopped 
at iteration n when \p["^ — p["~^^\ < e). It was found that decreasing e to 10"-^° did not change the results by more 
than 1% (but in turn the computational time was increased by a considerable factor). 

Consider now the variations Spi = pi{p + Sp) — pi{p) in the local fluid densities when increasing the chemical 
potential by a very small (but finite) increment Sp (for instance Sp = 10^^, with Wff taken as the energy unit). We 
show in Fig. 1 a small portion of a typical low-temperature isotherm computed in a 87% porosity sample of size L — 50 
at T* = T/Tc = 0.5 {ksTc/wff = 2). The random character of the condensation events results in the fine staircase 
structure of the curve, with the vertical steps corresponding to the avalanches separated by nearly horizontal portions 
corresponding to the smooth deformation of the grand-potential local minima (note the small scale on the vertical 
axis). The slope of these near-plateaus depends on temperature and its order of magnitude can be estimated from 
the isothermal compressibility of the pure lattice gas. It turns out that all the p^'s vary by a very small amount in 
these fiat portions of the isotherms (less than 10~^ — 10~^) whereas the vertical steps are due to a significant variation 
of some of the p^'s (a low temperature, this typically corresponds to a gazeous region becoming liquid). We can use 
this feature to identify the avalanches unambiguously. Specifically, we consider that the change Spf = J2i ^Pi in 
isotherm is only due to a smooth deformation of a grand-potential minimum when (5pmax = maxi{(5pi} is smaller that 
some threshold, which we take equal to 0.1. Otherwise there is an avalanche, and we associate to it all the sites for 
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FIG. 1: Portion of a typical adsorption isotherm in a 87% porosity aerogel at T* — T/Tc = 0.5 (p/ is the average fluid density 
inside the aerogel and jj. the gas chemical potential). Calculations were done in a sample of size L = 50 using the dichotomy 
procedure explained in the text. 



which 5pi > 0.001. The size of the avalanche is then s = J2i av ^Pi^ where the sum is restricted to these sites|2l| (by 
taking s = Spf = '^^Spi one would neglect the small increase of which is due to the "smooth" variation of the 
Pi's on the other sites of the lattice). This criterium thus implies that avalanches for which (5/3max < 0.1 cannot be 
detected. By choosing a smaller threshold, e.g. (Spmax < 0.01, one may find new avalanches, but their size s is very 
small (always less than 10). As a consequence, avalanches with s < 10 will not be considered in the following. In any 
experiment there is also a threshold below which signals cannot be detected. 

Since one must collect the results of many samples to get a good statistics of the avalanches properties, it is 
impossible in practice to study a whole isotherm with an increment 5p as small as 10~^. It is much more efficient 
to proceed by dichotomy, comparing the configurations obtained for two values of p, say pi and p2 (initially, pi is a 
large negative value and p2 = —4, the value at saturation). If 5/3,„ax > 0.1 between the two configurations, one then 
compares the configurations for pi and {pi + /X2)/2 on the one hand, and {pi + p-i)!^ and p2 on the other hand. 
One then looks for an avalanche in the two intervals, and the same procedure is repeated until no avalanche is found 
between two configurations, or the difference in chemical potential becomes smaller than 10~^ (this is the method 
that was used for obtaining the isotherm shown in Fig. 1, which explains that the points are irregularly spaced). 

Finally, to characterize geometrically an avalanche, which is a random object, one may define an effective "radius 
of gyration" Rg by attributing to each site i a "mass" that corresponds to the change in the local density. 



E»,aiX^V-r^)2(5pi 



Rl = ^^'"v, ;y " (3) 



where the sum is restricted to the sites contributing to the avalanche (i.e. with 5pi > 10 and r,; 
Si av ^i^Pi/^i av ^Pi "ccutcr of mass" of the avalanche E^. 



B. Number, size and shape 



Let us first study how the number of avalanches with a given size varies with temperature and porosity. To illustrate 
the general trend, we consider isotherms obtained in 87% and 95% porosity aerogels, as shown in Figs. 2(a) and 3(a). 
The sequence of avalanches was recorded in 500 samples of linear size L = 50 and 50 samples of linear size L = 100, 
respectively (this approximately corresponds to the same ratio L/^g « 10, where S.g is the correlation length of the 
aerogel l]). Figs. 2(b) and 3(b) show N{s), the (unnormalized) avalanche size distribution (/ N{s)ds — Nav, the total 
number of avalanches of size larger than 10). As T increases, Nav decreases from 335000 to 81000 and from 83700 to 
19800, respectively. The size of the largest possible avalanche corresponding to a simultaneous change of all the pi's 
from to 1 is 2L^. One thus has Smax — 2.5 10^ and 2. 10^ in the 87% and 95% aerogels, respectively. 

It is readily seen that there exist much larger avalanches in the lighter aerogel. This is a direct consequence of the 
larger void space available, as discussed in our previous papersP,^]. In particular, at T* = 0.5, some avalanches span 
the whole system, signaling that a macroscopic avalanche possibly occurs in the thermodynamic limit. Note however 
that the statistics for very large avalanches is not good, which explains the erratic behavior seen in Fig. 3(b) for large 
s. Moreover, finite-size effects play an important role in the large-porosity/low-temperature regime and they should 
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be taken into account to properly describe the evolution towards the thermodynamic limit (see Ref . jl^ for a careful 
study in the case of the RFIM). Unfortunately, this would require a considerable amount of numerical work, and we 
shall therefore only concentrate our discussion on the non-spanning avalanches. 




FIG. 2: (a) Average isotherms in a 87% porosity aerogel at T* = 0.5,0.65, and 0.8 (from top to bottom); data correspond to 
an average over 500 gel samples of size L — 50. (b) Corresponding (unnormalized) avalanche size distributions N{s). The solid 
lines represent the fit N{s) ~ s'^e"*''^" (Color on line). 




FIG. 3: Same as Fig. 2 in a 95% porosity aerogel at T* — 0.5, 0.8, and 0.9; data correspond to an average over 50 gel samples 
of size L = 100 (Color on line). 

Figs. 2(b) and 3(b) show that the evolution of N{s) with temperature is quite similar in the two aerogels: as T 
increases, both large and small avalanches are suppressed, and, as a consequence, the total number N^v decreases. 
Only the number of avalanches of intermediate size stays roughly constant with T. If one refers to what happens 
in the field-driven RFIM at T = 0, this result may appear surprising. Indeed, the isotherms get smoother as one 
increases the temperature, like the magnetization curves M{H) in the RFIM as one increases the width of the 
random field distribution |9l linl|. This seems to imply that the temperature plays the role of an effective disorder in 
the system. Therefore, one would expect that large avalanches become rare as T increases (which is indeed what 
happens), but also that small avalanches proliferate, so that Nav increases on average It turns out, however, 

that temperature also plays another role in our system, a role which is well documented in the literature devoted to 
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FIG. 4: Point cloud of the avalanche size s recorded in the 87% porosity samples as a function of chemical potential /i at 
T* = 0.5 (a) and T* = 0.8 (b) 

capillary condensation'23'|: adsorption in a small cavity becomes more gradual as T increases. Specifically, mean- field 
calculations for a single pore of simple geometry (for instance a slit of width I) predict that adsorption is continuous and 
reversible above some critical temperature Tc{l) which decreases as / decreases. Since small avalanches are associated 
with gas condensation in small cavities or crevices defined by the neighboring gel strands, this phenomenon explains 
the decrease of the number of small avalanches (when T > Th{(f>), the temperature at which the hysteresis disappears, 
there are no avalanches at all and the isotherms become reversibleQ)- It is clear that the random- field picture can 
only be applied at a scale much larger than the correlation length of the solid: schematically, capillary condensation 
prevails for s^^^ < whereas the disordered character of the aerogel structure plays a major role for s^^^ > ^g- The 
two mechanisms compete at the scale s^^^ ^ which explains that the number of avalanches of intermediate size 
does not vary significantly with T. 

If we now concentrate on the disorder-dominated regime, the most remarkable feature, quite visible in the 95% 
gel at T* = 0.5, is that N{s) follows an approximate power-law behavior over several decades. This is the defining 
signature of criticality, an issue that has been at the focus of intense research activity in recent years, for instance to 
explain the scaling properties of the Barkausen noise in magnetic materials!^. In the RFIM picture, the power-law 
behavior is related to the existence of a disorder-induced critical point 0,^3- this context, the avalanche size 
distribution should be analyzed using the form 

Nis) ^ s-"$(s/so) (4) 

where r is a critical exponent (25l|, $ is a scaling function, and sq a cut-off that depends on the distance to the critical 
point (sq diverges as the critical disorder is approached). Our results are however too limited to obtain sensible scaling 
collapses (although a crude attempt shows that the scaling function is non-monotonic, in agreement with the RFIM 
behavior (IQi]). We simply replace the function $ by an exponential form only describing the cut-off of the power-law 
at the scale sq. Figs. 2(b) and 3(b) show the fits of N{s) in the two aerogels according to this approximate treatment. 
The power-law regime is very limited in the 87% aerogel even at the lowest temperature, indicating that the system 
is clearly outside the critical region. On the other hand, in the 95% gel, the curves show a linear portion on a log-log 
scale that increases as T decreases: the cutoff sq is pushed to larger values of s (from 3200 to 29500 as T* varies 
from 0.9 to 0.8) and reaches the box size for T* = 0.5 (our set of data is however too limited to estimate the critical 
exponent associated to the divergence of sq). The exponent t increases and becomes approximately equal to 1.8 in 
the fully developed power-law regime at T* — 0.5. Taking into account the fact that our systems are small and that t 
in a finite system depends on disorder 19], this value does not appear inconsistent with that predicted for the RFIM, 
T = 2.03 ± 0.03J0J. In any case, the gcl-fiuid system is clearly close to criticality at this temperature. 

The number and the size of the avalanches along an adsorption isotherm strongly depend on the value of the 
chemical potential. We thus also recorded the value of fi at which each avalanche of size s occured. Figs. 4 and 5 
show the corresponding point clouds in the 87% and 95% gels at T* = 0.5 and 0.8, and T* = 0.5 and 0.9, respectively. 
As a general rule, the size of the largest avalanches increases with /i, especially in the last stage of the adsorption 
(note the logarithmic scale on the vertical axis). But the most striking feature is that avalanches of very different 
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FIG. 5: Same as Fig. 4 in the 95% porosity samples at T* = 0.5 (a) and T* = 0.9 (b). 



sizes occur in a small interval of /i. For instance, in the 87% gel at T* = 0.5 and ~ —4.5, s varies from 10 to 
more than lO'^. In the 95% gel at the same temperature and fi « —4.12, the avalanches span almost all possible sizes 
between 10 and 10^. These observations are in complete contradiction with the picture based on the independent pore 
model which is often used to interpret capillary condensation Dhenomenaj2fi|| . In this model, only pores of a certain 
characterisic size I fill at a given value of the chemical potential, and the adsorbed fluid density is directly related 
to the number of pores of this size that can be found in the material. (One then try to extract from the isotherm 
the pore size distribution.) This picture is clearly wrong in the present case because of the very open and intricate 
character of the gel structure that allows condensation events of very different sizes. 

The geometrical characterization of the avalanches is rather difficult. Using s^^^ as an estimate of the extension of 
an avalanche is misleading because events with very different shapes are observed (see for instance Figs. 6 and 7 in 
Ref.0). This variety in the geometry is reflected in the values of the radius of gyration as defined by Eq. (3). 
The value of Rg is indeed proportional to the maximal extension of the avalanche, but the prefactor depends on its 
shape. Figs. 6 and 7 show how Rg varies with s in the 87% and 95% aerogels. One can see that avalanches with a 
similar size may have quite different radii of gyration. For instance, in the 87% aerogel at T* = 0.5, avalanches of 
size s ^ 10'^ have a radius that varies between 4 and 10 (note that for a compact object, a factor of 2 in Rg would 
correspond to a factor 8 in s). Similarly, in the 95% aerogel at the same temperature, Rg varies from 10 to 15 for 
avalanches of size s ~ 10*. It is therefore impossible to associate a unique characteristic length to an avalanche of a 
given size. Moreover, by examining the avalanches that occur in a small interval of /i, we find that they also have 
very different radii of gyration. 

As the temperature increases, the values of Rg for a given size s are less scattered. Indeed, avalanches become 
more compact as they correspond to condensation events that occur in a single cavity of the aerogel. This evolution 
is very clear in the 87% aerogel: at T* — 0.8, the radius of gyration of most avalanches is between 2 and 4 and there 
are almost no avalanches with Rg larger than 6. On the other hand, there is still a large domain of variation of Rg 
in the 95% aerogel at T* = 0.9, and one should go to even higher temperatures to observe a significant reduction. 
There is in fact a remarkable resemblance between Figs. 6(a) and 7(b) which suggests that the avalanche properties 
are similar apart from a scale factor of about 2. The adsorption isotherms have also similar shapes, as can be seen in 
Figs. 2(a) and 3(a). This confirms the observation made in Ref.^f that the adsorption processes in a low porosity gel 
at low temperature and in high porosity gel at a higher temperature have similar properties. 

C. Spherical avalanches 

In both Fig. 6(a) and Fig. 7(b), one may notice that a small set of avalanches stands out . In the 87% gel, at 
T* — 0.5, these avalanches have a size s in the range 100 — 500 and a smaller radius of gyration than the other 
ones. They are therefore more compact, and a direct observation of some of them shows that their shape is almost 
perfectly spherical. This is confirmed by a systematic study of the asphericity parameter which is always smaller 
than 10~^i27|. As shown in the inset of Fig. 6(a) on a log-log scale, the radius of gyration of these avalanches 
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FIG. 6: Radius of gyration Rg and size s of the avalanches in the 87% porosity samples at T* = 0.5 (a) and T* — 0.8 (b). The 
inset in (a) is a magnification of the point cloud of the spherical avalanches (see text) in log-log scale; the solid line represents 
the equation Rg = 0.42 s'''^^ (Color on line). 




FIG. 7: Same as Fig. 6 for the 95% porosity samples at T* = 0.5 and 0.9. In the inset of (b) the solid line represents the 
equation Rg = 0.73 s"'^^ (Color on line). 



varies approximately like Rg « 0.42 s" '^^ which is indeed close to the expected value for spherical compact objects 
(Rg = ^/3J5 R = [3s/(47r)]i/3 0.485^/3 for a sphere of radius R). The radius R varies between 2.5 and 4, this 

latter value corresponding approximately to which is also more or less the size of the largest voids in the aerogel 
Moreover, as shown in Fig. 8, most of these avalanches occur at the end of the adsorption process: they thus 
correspond to the condensation of the last gas bubbles that remain in the large cavities. 

A similar behavior is observed in the 95% aerogel at T* = 0.9 (Fig. 7(b)). The avalanches are again approximately 
spherical and occur in the last stage of the adsorption process. However, since the temperature is higher, the boundaries 
of the gas bubbles are less well defined, and the radius of gyration of varies approximately like Rg « 0.73 s°-^^. The 
largest observed value of i? = ^5/3 Rg is again comparable to the size of the largest cavities in the aerogel (^g ^ 10) . 
Although there are also some other remaining bubbles of more complicated shape (see Figs. 7 and 9 in Ref.0]), 
these spherical avalanches actually represent the only case where one can relate directly the characteristic size of the 
condensation events to the structure of the aerogel. We stress again that, in general, the avalanches have a complicated 
shape so that one cannot describe the condensation of the adsorbed gas in terms of "constant curvature bubbles" 0| . 
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FIG. 8: Probability distribution P{fi) of spherical avalanches along the adsorption isotherm in the 87% porosity aerogels at 
T* — 0.5. The average isotherm is shown as a continuous line. 




FIG. 9: (a) The undersaturation A/i — fisat — pi as a function of the inverse radius of the spherical avalanches in the 87% and 
95% aerogels at T* — 0.5 and T* = 0.9, respectively, (b) is a magnification of the small point cloud in (a) near the origin 
corresponding to the avalanches in the 95% aerogel. The solid line with a slope 1 corresponds to the Kelvin equation (see text) 
(Color on line). 



The fact that spherical avalanches are clearly identified gives us the opportunity for testing a macroscopic description 
of the gas condensation in terms of a competition between volume and surface contributions. Such a description 
underlies the classical interpretation of capillary condensation based on the Kelvin eauation|26l| that is commonly 
used to extract from the isotherms information on the characteristic pore size in the solid. We thus write the change 
in t he g rand potential associated to the condensation of a spherical gas bubble of radius R at the chemical potential 
fi as [3^ 

A^^i^l)^[pi^i)-p+i^l)]—R'-JinR^ (5) 

where is the corresponding bulk pressure, p^{fJ,) is the pressure of the metastable liquid at the same chemical 
potential (both pressures are computed from the equation of state of the bulk lattice gas), and 7 is the (planar) 
liquid-gas surface tension. All these terms except R are temperature dependent. Assuming that condensation occurs 
at equilibrium, the Kelvin equation (in its Laplace form) corresponds to setting AQ = in the equation. Expanding 
as usual p and about Psat and to first-order in A/i = fisat — (which is valid for small undersaturations) , this yields 
A/j = Sj/Kpi — pg)R\, where pi and pg are the densities of the coexisting liquid and gas at this temperature. Since the 
surface tension on a lattice depends on the orientation of the surface with respect to the lattice vectors, we estimate 
7 by taking the average of the three values calculated in the principle lattice planes (100), (110) and (111)0|. The 
corresponding "Kelvin plot" is shown in Fig. 9. We see that the "macroscopic" Kelvin equation provides the correct 
trend for the variation of A/i as a function of in spite of the fact that actual condensation does not occur at 
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equilibrium {Ail ^ 0(23) and that the adsorbent potential and the curvature corrections are neglected in Eq. (4) 
(note that is ten times smaller in the lighter gel, which is due to the effect of temperature, correctly described by 
the term 2>^/{pi — Pg), and to the change in the average radius i? by a factor of about 2). The scatter in the data, 
however, shows that the Kelvin equation is only an approximation. 

IV. SCALING OF THE ADSORPTION ISOTHERMS AND PHASE DIAGRAM 

The existence of a critical point in our system, at which an infinite avalanche occurs for the first time in the thermo- 
dynamic limit and the avalanche size distribution has a power-law behavior, suggests the possibility of performing also 
a scaling collapse of the adsorption isotherms. In the case of the nonequilibrium RFIM at T = where the disorder is 
controlled by the width a of the distribution of the random fields, standard renormalization group arguments suggest 
that the magnetization curve M (H) for different values of a may be collapsed using the scaling form(l(]|| 

M{<j, H) - M,(a„ H,) ^1 r f M±{h/ \ r f') (6) 

where Mc is the critical magnetization (i.e. the magnetization at the critical field He for the critical disorder CTc), 
r — {<7c — o')/ac and h — H — He are the reduced disorder and reduced field respectively and M.± is a universal 
scaling function (± refers to the sign o f r) . The estimated values of the critical exponents in three dimensions are 
P = 0.035 ± 0.028 and pS ~ 1.81 ± 0.320. 




FIG. 10: Average isotherms in aerogels of varying porosity (the fluid density p/ is rescaled by the density at /i = fisat)- From 
left to right: 4> = 87%, 90%, 92%, 95% at T* = 0.4, 0.5, 0.6, 0.7, 0.8. The curves result from an average over 500 samples of size 
L = 100 (Color on line). 

The situation in our system is more complicated than in the RFIM because the Hamiltonian (or, more correctly, 
the grand-potential functional fl) contains three parameters that can be tuned independently (the porosity 0, the 
interaction ratio y, and the chemical potential fi), and, moreover, the temperature controls the ruggedness of the 
free-energy landscape. Although a rigorous proof is lacking, it seems however reasonable to assume that the critical 
behavior is governed by a unique fixed point, which may or may not be the same as for the T = nonequilibrium 
RFIM. In other words, there are only two relevant variables from a renormalization-group standpoint: the chemical 
potential (that plays the role of the magnetic field H in the RFIM) and an effective disorder a which is some smooth 
funtion of y, (j) and T. (In addition, one could also have to take into account the microstructure of the aerogel, on top 
of the dependence on the porosity: in particular, one may suspect that the fractal character of very light aerogels, 
e.g. (j) > 99%, modifies the values of the critical exponents or, at least, induces a different scaling behavior in some 
intermediate regime; this possibility is not taken into account in the present study which does not consider aerogels 
of porosity larger then 95%.) For fixed y and T, there is a critical value <^c(?/, T) of the porosity, and a critical value 
pLc{y,T) of the chemical potential. Alternatively, for fixed y and cj), there is a temperature Tc{y,<j)) and a chemical 
potential /j,c(2/, ^) at which the system becomes critical. If this picture is correct, one expects to observe the same 
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critical behavior whether the critical line 4>c{T) (or Tc{(j))) is approached along a constant-(/) or a constant-T trajectory 
in the temperature-porosity diagram (for simplicity, we forget the dependence on y which is kept constant and equal 
to 2 in all our computations 

Consider for instance the collection of isotherms displayed in Fig. 10 which covers a significant range of temperatures 
and porosities. This would correspond to 2.08K < T < 4. 16 A' for the real ^He-aerogel system, but this must be 
considered only as a crude estimation: we indeed recall that our treatment is not aimed at a quantitative description of 
actual systems because of the coarse-grained character of the model and the use of a (local) mean-field approximation. 
Note that the density of the adsorbed fluid is rescaled by its value a,t fi = figat so that all the curves go to 1 as yit — > fisat 
(one has pf{(j),^sat) ~ (f' Pi{T) where pi(T) is the density of the bulk liquid at coexistence). 
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FIG. 11: Scaling plot of the slope of the T* = 0.4 isotherms for (j) = 87%, 90% and 92% according to Eq. 8 where x = 
[M-Mma.]/ I (l>-MT) and y =| - 0c(r) f"-'^ pf{(l), Hsat)-^dpf{(j,, fi) /dfi with /3 = 0.1, I3S = 1.7, and ,^^(0.4) = 0.943 
(Color on line). 

Let us first study the scaling behavior when varying the porosity at fixed temperature. By analogy with the RFIM, 
our basic scaling hypothesis is the existence, in the limit L oo, of only two scaling variables u{4), p) and v{(j3, p) 
that measure the distance to the critical point {(j)c{T) , pdT)) . Accordingly, the adsorption isotherms in the vicinity 
of the critical point {(j)c{T), Pc{T)) should scale as 

Pf{(i),p) - pf{(i)c{T),Pc{T)) - pf((t>,p,at) I u/uo f j/j^^ps ) C^) 

where TZ± are the corresponding scaling functions and uq and vq are two nonuniversal scale factors that are 
temperature-dependent (they actually depend on all the specific features of the model js^]: y, the lattice structure, 
the range of interaction, etc.). Note that the scaling variables u and v need not be exactly (/> — (j)c{T) and p — pc(T), 
but their dependence on (j) and p should be analyticj^. Unless one uses the proper variables, the curves do not 
collapse. In the case of the RFIMT^, it was shown that the proper scaling must be done with replacing h — H — He 
in Eq. (6) by the "rotated" variable h' — h + br, where b is some non-universal constant (on the other hand, there 
is no need to rotate r because the corresponding correction is irrelevant when a —> Uc)- In practice, the rotation of 
the X-axis can be avoided by simply lining up the peaks in the derivative of the curves, with no correction terms, as 
is often done in analyzing experimental data (by considering the slope dpf((j),p)/dp instead of pf{4>,p) we also get 
rid of the dependence on pf{(j)c{T),pc{T))). What to do on the y-axis is a matter of choice and, for simplicity, we 
just take u = (j) - (pdT) . (Note that in Ref.fT^ the curves dM/dH in the RFIM are plotted using r — [a — (Jc)/cr 
instead of the more standard choice r ~ {a — ac)/o'c', alternatively, it has been proposed to expand r to second-order 
in (a — a^) /a^\l^: these different choices are equivalent asymptotically if no rotation is required, but they may change 
significantly the amplitude of the critical regions.) Therefore, to summarize, we look for a scaling of the form 




(8) 
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where u = <j) — (f)c{T), /Xmaa; T) is the value of the chemical potential at the maximum of dpf{(f>, 'R,± is the 

derivative of Ti±, and a(T) and b{T) account for the dependence on uq and Vq. Since it is important to have smooth 
isotherms so that one can take the derivative numerically, the average must be performed over a significant number 
of samples: the curves shown in Fig. 10 result from an average over 500 gel realizations of size L = 100 (however, 
there are still some annoying irregularities in the average slopes, especially around /i = fJ-max)- 

Calculations are of course done in finite systems so that one should also take into account finite-size effects 0. These 
are especially important in the low-disorder regime {((> > 4>c{T)) in which there is a discontinuity in the isotherms in 
the thermodynamic limit. This introduces an additional difficulty in the analysis and it is a very demanding numerical 
task to perform a systematic finite-size scaling study of all the curves (see Ref.0 for such a study of the T* = 0.45 
isotherm in the 95% porosity aerogel). As a consequence, we shall only consider here the isotherms that are clearly 
in the strong disorder regime, which leads to reasonably size-independent results: this eliminates in particular the 
T* = 0.4, 0.5 and 0.6 isotherms in the 95% porosity aerogel. 

Fig. 11 shows the scaling plot of the slope of the T* = 0.4 isotherms for (j) = 87%, 90% and 92% obtained with the 
choice (36 = 1.7, /3 = 0.1, and ^c(0.4) = 0.943 (the scale factors a{T* = 0.4) and b{T* = 0.4) are arbitrarily fixed 
equal to 1). It is actually very difficult to adjust independently the three unknown parameters (36, (3 and (f>c{T); the 
values of the exponents used here are the ones that allow for a collapse of a whole set of isotherms, as discussed below. 
It can be seen that the (/> = 92% and — 90% isotherms collapse onto a single curve, which is not the case of the 
(f> = 87% isotherm, especially on the left-hand side of the peak. This part of the isotherm is indeed dominated by the 
presence of a liquid layer that coats the aerogel strands PjQ, and in the low-porosity gels the remaining available void 
space is too small for allowing large avalanches to spread out. This contribution to the adsorption has nothing to do 
with a critical phenomenon and it is reasonable to eliminate from the scaling plot all the isotherms obtained in the 
87% porosity aerogel. More generally, at each temperature, we choose to discard the isotherms that do not belong 
to the scaling region because the porosity is obviously too far from 4>c{T) (i.e., the distance from the effective critical 
disorder is too large). This finally leads to the global scaling plot shown in Fig. 12, obtained by adjusting also the 
nonuniversal scale factors a(T) and h{T) (which, however, always remain of order 1). 

The overall quality of the collapse is not perfect, but similar to that found in the RFImT^ (it is actually better 
because we only keep two isotherms at each temperature 32j). In fact, as discussed in Rcf. 10], the magnetization 
curves (or, here, the isotherms) are quantities that do not collapse well, which makes it difficult to obtain from the 
scaling plots a precise determination of the critical exponents or of the location of the critical point. Although our 
estimated values of the exponents (35 and (3 are compatible with those of the RFIM, the uncertainty is much too large 
to conclude that our model belongs to the same universality class. On the other hand, it is noteworthy that the value 
1.7 for the exponent (36 is consistent with the value 1.22 for the finite-size scaling exponent x = (^^/^ obtained in 
Ref.0, if one uses for l/i^ the RFIM value, 0.71 ± 0.09[l3. 
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FIG. 12: Scaling plot of the slope of the isotherms according to Eq. 8 (when varying the gel porosity at different temperatures) 
or to Eq. 9 (when varying the temperature for the porosity (f> = 95% ). x and y are the corresponding scaling variables and 
scaling functions (Color on line). 
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FIG. 13: Approximate phase diagram of the model in the temperature-porosity plane deduced from the scaling plots. The 
dashed curve is the critical line Tc{(j>) (or <j)c{T)) below which there is a first-order nonequilibrium transition in the adsorption 
isotherm. The squares indicate the isotherms that are included in the scaling plot of Fig. 12 (Color on line). 
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FIG. 14; Scaling plot of the slope of the T* = 0.8 and T* = 0.7 isotherms in the 95% aerogel according to Eq. 9 where 
X ^ ii^aAl I r/rc(0) - 1 I"* and y =| T/re(0) - 1 I"*-" p[{(l},ti,aty^dpf{(l,,ti)/d^ with P = 0.1, (55 = 1.7, and 

T*{(f) = 0.95) = 0.5 (Color on line). 

From the values of 4>ciT) extracted from the scaling plots at each temperature, we can map out an approximate 
phase diagram of the model in the temperature-porosity plane. This is shown in Fig. 13, where the critical line 
has been forced to extrapolate towards the zero-disorder limit <f>c{Tc) = 1- (Let us recall however that, in a real 
system, there cannot be a clear-cut transition close to Tc because thermal fluctuations become larger and can allow 
the system to cross the free-energy barriers on the experimental time-scale; moreover, as noted before, we expect the 
fractal character of the microstructure in very light aerogels to come into play.) The present calculations indicate 
that the critical temperature drops very quickly with the porosity, and the extrapolation to T = predicts that there 
is no transition at finite temperature in gels of porosity smaller than (pciT = 0) « 92.5%. This is in agreement with 
finite-size scaling calculations performed in the 87% porosity aerogel that show that the isotherms are continuous at all 
temperatures in the thermodynamic limit. On the other hand, in the 95% porosity aerogel, the study of Ref.Q predicts 
a transition at T* — 0.45, which is indeed below the critical line of Fig. 13. Morever, as shown in the preceding section, 
the avalanche size-distribution has a power-law behavior at T* = 0.5, which is the critical temperature extracted from 
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the scaling plot {T*{4> = 95%) = 0.5 ± 0.02). It turns out that the same critical temperature is also obtained by 
studying the scaling behavior when varying the temperature at fixed porosity. The scaling is now performed according 
to the expression 

— a{(f>)\u\P n±[b{(f>) j (9) 



where u = {T — Tc{(l)))/Tc{(j>) and a'(0) and 6'(0) are porosity-dependent scale factors. It is clear from Fig. 13 that this 
procedure can only be used with the 95% porosity aerogel. The scaling plot of the T* = 0.8 and T* = 0.7 isotherms is 
shown in Fig. 14 (again, the T* = 0.6 isotherm is discarded because of the presence of important finite-size effects). It 
can be seen that a good scaling is obtained with the same values of the critical exponents as used before. Remarkably, 
the resulting scafing function is also the same and it can be included in Fig. 12 by adjusting the scale factors a'{4>) 
and b'{(p). This is a clear indication that the critical behavior of the system is controlled by the same fixed point for 
all temperatures and porosities. 

We now discuss the comparison of all these results with experiments. 



V. DISCUSSION 

In this work and our earlier workQ, Q, we have presented a detailed numerical study of the adsorption of "^He 
in silica aerogel that explains the modification in the shape of the isotherms as the temperature or the porosity is 
changed. We have found that the adsorption occurs via a sequence of irreversible condensation events (avalanches) 
of varying size and shape. The traditional picture of capillary condensation in mesoporous solids in which pores of 
increasing size fill up successively as the chemical potential or the gas pressure increases does not hold in aerogel, 
except, approximately, for the special category of spherical avalanches that occur at the end of the adsorption process. 

Our calculations predict that there exists a line of critical points in the temperature-porosity diagram separating 
two regimes of continuous and discontinuous adsorption and characterized by a power-law behavior of the avalanche 
size distribution. There is a rather large scaling region where adsorption isotherms at different temperatures and 
different porosities collapse onto a single universal curve. Although the system sizes considered in this work are too 
small to reach a definite conclusion, the estimated values of the critical exponents seem to be consistent with those 
of the three-dimensional field-driven RFIM at T = 10, 19]. 

The qualitative agreement between the simulated and experimental loops is a strong argument in favor of the 

validity of the above scenario. However, a complete validation still awaits the experimental description of the scaling 
behavior of the isotherms and the observation of a power-law behavior in the avalanche size distribution. This latter 
issue requires some additional comments. It may be indeed distressing that avalanches have not been observed so far 
with *He in aerogel. In fact, the only exemple of avalanche behavior in adsorption/desorption experiments concerns 
the draining of superfluid ''He from Nucleporej^Sli a porous material that contains a random spatial distribution 
of interconnected cylindrical pores. In this case, superfiuidity appears to be a critical factor for the existence and 
propagation of avalanches via the presence of a superfiuid film on the Nuclepore substrate that provides a coupling 
mechanism even between distant pores. It is likely however that such a mechanism is not necessary for avalanches to 
occur in aerogels because of very open character of the microstructure. 

There are actually several requirements for observing avalanches in sorption experiments, which are all satisfied in 
the case of the experiments with Nuclepore 34] , but not in the case of the existing aerogel studies 0, |^ IE IS- Let 
us first recall that, by definition, avalanches correspond to some collective reorganization of the system in response 
to a variation of the external driving force. If one is changing the density rather than the chemical potential or the 
pressure (by introducing a controlled amount of gas into the sample and then measuring the pressure change), the 
system may follow a different path among the metastable states, and the vertical jumps in the isotherms that are the 
signature of the avalanches may not be observed 35] (this could explain that there are many data points along the 
adsorption isotherm at T = 2.3AK in 98% aerogel displayed in Ref. 4]). Secondly, as already noted, the observation 
time must be larger than the time the system takes to respond to the change in the external control variable. This can 
be done by continuously changing the chemical potential at a very slow rate, and measuring the amount of gas in the 
solid regularly (every x seconds). In practice, this procedure is only possible if the (local) equilibration time is short 
(for instance, it is smaller than one second with superfluid helium in Nuclepore). Unfortunately, the relaxation times 
in aerogel are much longer (of the order of two hours with '*He at low temperature in the experiments of Ref.0), and 
one has to change /i or P by a finite step and keep it fixed until a new local equilibrium (i.e., a new metastable state) 
is reached. With this procedure, the issue of having a system with a long equilibration time boils down to stability of 
the experimental conditions and of course patience; in principle, one could then observe avalanches. Yet, if A/i or AP 
is not small enough, the avalanches may overlap and the resulting isotherm may appear continuous. This is probably 
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the situation with the experiments in Ref.|^ where AP/P ~ 10~^ (to be compared with AP/P = A^/{kBT) ~ 10~^ 
in the Nuclepore experiments [s^])- Finally, a last requirement is that the sensitivity of the measurements is large 
enough to discriminate the jumps in the adsorbed fluid density. This condition is more easily satisfied when these 
jumps are large enough, which implies that one must be close to the critical line displayed in Fig. 13. Although one 
must be very cautious in comparing our temperature scale with the true experimental one, it seems that this last 
condition is also not satisfied in the existing experiments. 
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